Real-world options: smile and residual risk 



Jean-Philippe Bouchaud 1 , Gulia Iori 2 and Didier Sornette 3 
1) Service de Physique de I'Etat Condense, CEA-Saclay, 91191 Gif sur 

Yvette CEDEX, France 

2) Service de Physique Theorique, CEA-Saclay, 91191 Gif sur Yvette 

CEDEX, France 

3) Laboratoire de Physique de la matiere Condensee, Universite de 
Nice-Sophia Antipolis, B.P. 70, Pare Valrose, 06108 Nice CEDEX 2, 

France. 



Abstract 



We present a theory of option pricing and hedging, designed to address 
non-perfect arbitrage, market friction and the presence of 'fat' tails. An im- 
plied volatility 'smile' is predicted. We give precise estimates of the residual 
risk associated with optimal (but imperfect) hedging. 
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1 Introduction 



In 1973, Black and Scholes []IJ developed the mathematical theory of perfectly 
hedged options in an arbitrage-free lognormal random walk asset price model. 
Fundamental for the mathematical formulation and deeply rooted in the mind 
of both theoreticians and practitioners, evolved the notion that a general 
strategy of pricing and coverage of derivatives involves the hedging away of 
all risks. However, we do not live in a Black-Scholes world : 

• asset prices do not follow a continuous time lognormal processes, but 
may exhibit large jumps and be distributed with fat tails, in addition to 
present subtle correlations between increments at different times (the 
markets are 'incomplete'), 

• the replicating portfolio cannot be implemented exactly, since it in- 
volves continuous rebalancing: several 'imperfections' such as trans- 
action costs, delays and lack of liquidity necessarily introduce a finite 
time scale for transactions. 

Confronted with these practical problems, many authors have attempted 
to generalize the Black and Scholes strategy. Here, we present an intuitively 
appealing formalism which is flexible enough to address these problems ef- 
ficiently and show how it can be implemented practically. Our method is 
based on a global (integral) representation of wealth balance, in contrast to 
Black-Scholes local (differential) approach. The mathematical tool-box is 
that of functional integration and derivation, in contrast to the standard Ito 
calculus and partial differential equation formalism. We propose two main 
new concepts : 

• The minimization of the non zero residual risk as a criterion to fix an 
optimal strategy [0. One should note in this respect that somewhat 
related ideas have previously been proposed in the mathematical liter- 
ature in a rather formal way . Furthermore, our hypothesis are more 
general and the resulting formulae are easily implemented numerically. 

• For strongly fluctuating assets, option pricing and its hedging is ob- 
tained by a 'tail chiseling' of wealth distributions as a treatment of 
large risks [§]. 
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One of our main thrust is to demonstrate how these ideas can be put in 
practical terms, using relatively simple mathematics. 

Let us discuss it for the basic European call option pricing problem - 
generalisation to more complicated products does not involve any conceptual 
difficulties. Suppose that an operator wants to buy a given share, a certain 
time t = T from now (t — at which the current value is x(t — 0) = Xq), at 
a fixed 'striking' price x c . If the share value at t — T, x(T), exceeds x c , the 
operator 'exercises' his option, with an immediate profit difference x(T) —x c . 
On the contrary, if x(T) < x c the operator may not buy the share. What 
is the price C(xq,x c ,T) of this possibility given to the operator by -say- the 
"bank" , and what trading strategy should be followed by the bank between 
now and T, depending on what the share value x(t) actually does between 
t — and t = Tl These are the two questions answered by Black and Scholes 
in the context of an ideal complete market, and which we wish to consider 
from a somewhat different point of view adapted to real-world markets. 

In common with Black and Scholes, the starting point is to write a wealth 
balance for the bank. However, instead of computing the instantaneous vari- 
ation of the value of its portfolio (which assumes continuous time), we write 
a global balance at time T, and assume that time is discrete, with a cer- 
tain 'microscopic' time scale r below which trading is impossible (or, in the 
presence of transaction costs, very unfavorable - see below). 

We shall assume for simplicity that the risk- free interest rate r is constant; 
more involved developments (including a random interest rate) are reported 
elsewhere. The global balance AW for the bank between t — and t = T is 
thus given by: 

AW = C(xo, x c , T) exp(rT) — max(x(T) — x c , 0) + 

T 

+ 4>(x,t)exp(r(T -t))[A T x -rx(t)] (I), 
t=o 

with x c the striking price, Xq the initial (t = 0) share price, and <f)(x, t) is 
the number of shares owned by the bank at time t, knowing that the share 
price is x. A T x is the difference of share values between times t + r and 
t. C(xo,x c ,T) is the looked for option price. Expression (1), which in the 
continuous case is in fact the time-integrated version of Black and Scholes's 
differential equation, has a very intuitive meaning : 
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• The first term is the gain from pocketing from the buyer the option 
price at t = 0, discounted at time T. 

• The second term gives the potential loss equal to — (x(T) — x c ) if x(T) > 
x c (i.e. if the option is exercized) and zero otherwise. 

• The third term quantifies the effect of the trading between t — and 
t = T: the true variation of wealth W between t and t + r is only due to 
the fluctuations of the share price, i.e. (f>(x, t)A T x, corrected by the fact 
that xtp(x, t) has not benefited from the risk-free interest rate. (Note 
that the term xA T (f> describes conversion of shares into other assets or 
the reverse, but not a real change of wealth). 

The bank wealth variation AW depends a priori on the specific realization 
{x(t)}\ t =o-T of the asset price. For a continuous log-normal process, the 
result of Black and Scholes is that AW is strictly vanishing for r — > for 
a particular choice of <f)(x,t) and C(xq,x c ,T). Equivalently, both equalities 
(AW) = and ((AW) 2 ) = hold, where (...) denotes the average over 
all possible different realisations of the history {x(t)}. In other words, the 
averages (...) are taken over the initial historical measure in which the drift 
is not in general equal to the risk-free interest rate. The first condition 
(AW) = is the 'no free lunch' condition. The second expression simply 
writes that the square of the bank wealth volatility (AW 2 ) is zero which 
ensures the hedging away of all risk. This is just another view point to 
retrieve Black and Scholes results, holding true for any quasi-Gaussian, time- 
continuous processes @. 



2 Risk-corrected option prices 

In the real world where AW cannot be made to vanish exactly, we propose 
to apply the spirit of Markowitz portfolio approach to the total bank wealth 
balance, viewed as an effective portfolio estimated at the time T of the death 
of the option. In this goal, one has to calculate the average return (AW) 
and the second moment (related to the square of the bank wealth volatility) 
(AW 2 ) when they exist (the case of very strongly fluctuating 'Levy processes' 
with ill-defined theoretical volatility will be addressed below). For the time 
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being, we assume that both (AW) and (AW 2 ) exist. The optimal strategy 
<p*(x,t) that the bank must follow should certainly be such that the uncer- 
tainty on the outcome is minimum. The reason for this is that the bank will 
adjust the option price according to its risk aversion, for instance with the 
criterion 

(AW)\ <M r = \VK : , (2) 

where 1Z* = (AW 2 ) — (AW) 2 |</,=0* is the square of the volatility of the bank 
wealth. Eq. (2) simply means that the bank wishes to keep the probability 
of global loss incurred by delivering the option at a certain level - say 10%. 
The minimisation of 1Z* is thus necessary to keep the option price itself as 
low as possible. Alternatively, the knowledge of (AW)]^^* and y/TZ* allows 
one to extract from the market price of a given option the value of its 'A' (a 
yet unused Greek letter). A can be thought of as an objective (dimensionless) 
measure of the real option price, in units of the real risk associated to the 
underlying asset - high A's corresponding to expensive options. Note that 
the term in the right hand side of eq.(2) allows the bank to introduce a 
bid-ask offer range around a central 'fair' price which itself is independent 
of A. This bid-ask offer range is dependent upon the bank portfolio and its 
risk-aversion. 



2.1 Independent increments of the share value 

Let us suppose that the local slopes A T x are statistically independent (but 
not necessarily Gaussian) for different times. In this case, both (AW) and 
the 'risk' measured as K[<f>(x,t)] = ((AW) 2 ) - (AW) 2 can be explicitely 
calculated. We can then determine the optimal strategy <p*(x,t) through a 
'functional minimisation': 

from which a general expression of both </>*(x,t) and of the residual risk 
1Z* = 7Z[4>* (x , t)} can be derived || - see the 'Technical sheet' below). These 
results are valid for an arbitrary stochastic process with uncorrelated incre- 
ments, including 'jump', or discrete-time, processes. The process can fur- 
thermore be explicitely time-dependent, with a variance which is a function 
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of time, as for the much studied 'ARCH' processes, which are Gaussian pro- 
cesses with a t-dependent variance V{t). 

The obtained formulae can be simplified in several cases [T2] and allows us 
to retrieve Black and Scholes' result when A T x are Gaussian variables and r 
tends to zero (continuous process). In particular, one finds indeed that TZ* 
vanishes exactly in the continuous time limit. One may also check that in 
that case, the results on the option price and the optimal strategy are indeed 
independent of the average return of the share. 

The condition TZ* = is very specific to the Gaussian can and does 
not hold for a more general stochastic process. This is the main difference 
between the present approach and that of Black and Scholes and subsequent 
workers (see however [[J): 

1) we find that a vanishing residual risk cannot be achieved in the general 
case and are able to quantify it precisely (see Fig 2 below); 

2) however, this does not imply that an optimal strategy does not exist. 
We have indeed found an optimal <p*(x, t) which minimize the risk and which 
is a natural generalization of Black and Scholes result. 

Note finally that other definitions of the risk are possible, for example 
through higher moments of the distribution of AW. The functional min- 
imisation technique presented here can easily be adapted to these cases - 
although the calculations are more cumbersome. 

2.2 Illustration: options on the MATIF 

These findings are relevant to various concrete situations. In particular, 
strong deviations from a Gaussian behaviour ( "leptokurtosis" ) are often ob- 
served in many situations. In this case, the above method allows one to 
estimate quantitatively the residual risk and correct the trading strategy 
and the option price accordingly. We present an illustration by compar- 
ing the option price obtained from the proposed functional quadratic mini- 
mization procedure and from Black and Scholes formula applied on MATIF 
options, with a 30 days maturity (T = Jjjr days). In this goal, we have deter- 
mined the historical conditional distribution P(x' ,t'\x,t) from the MATIF 
daily quotation variations in the period 1990-1992. We have tested a certain 
degree of stationarity of the data by constructing the distribution of price 
increments over various time intervals. We then determined P(x' .t'\x,i) by 
scanning the data set and using the assumption of stationarity implying that 
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P(x', t'\x, t) = P(x' — x,f — t). In order to apply the Black and Scholes re- 
sults, we have fitted P(x' — x, t' — t = 1) by a log-normal distribution, finding 
a volatility a 2 ~ per day. P(x' — x,t' — t = 1) and its log-normal fit 
is shown in Fig.l. We can also apply directly our formulae using the empir- 
ically determined distribution and obtain the option price C(xq, x c ,T) as a 
function of the striking price x c (Fig 2) in a risk insensitive world (A = 0), 
and compare it to the Black-Scholes price with the historical volatility (dot- 
ted curve). As is well known, this procedure underestimates the 'true' price: 
the presence of 'tails' in the distribution induces a larger effective volatility. 
One can in fact invert the Black-Scholes formula and determine an implied 
volatility from our determination of the price. The result is given in the in- 
set of Fig 2 and reproduces the well known 'smile': the effective volatility is 
stronger for 'out of the money' options. More importantly, the residual risk 
1Z* is non-zero, it furthermore depends on the trading 'frequency', i.e. the 
number of 'inactive' days r without rehedging. As expected, the larger this 
number of days, the larger the risk : we have plotted in Fig 3 the quantity 



ylZ* as a function of r, together with the total transaction costs K, associ- 

ated with these rehedging, assuming that each trading has a cost equal to 
0.05% of the share's value. The full curve represents the sum of these two 
sources of extra-cost: /C = + /C is the risk- aversion cost (A = 1) plus 
the transaction costs. Interestingly, K has a minimum, which is the optimal 
trading rate (r ~ 10 days for this particular choice of parameters). Due to 
the impact of transaction costs, it appears that a daily rehedging strategy is 
not reasonable - unless A is very large (strong risk aversion). Note that /C 
represents an appreciable fraction of the 'fair' price. 

Finally, we have compared in Fig 3 (Inset) the residual risk for our optimal 
strategy with the residual risk obtained following the Black-Scholes strategy. 
As expected, the latter is larger: the risk can be substantially reduced using 
our 'optimal' strategy. 

2.3 Rare events and fat tails 

Our method can be adapted to many situations: more exotic options, or 
more complicated stochastic processes, such as correlated Brownian motions 
[0] or Levy processes, which have been argued by many authors || to be 
adequate models for short enough time lags, when the kurtosis is large. The 
latter case is interesting since the notion of variance becomes ill-defined due 
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to the presence of extremely large fluctuation (crashes). The above criterion 
fixing the optimal strategy (f>(x,t), based on a minimization of the variance 
is thus meaningless. Intuitively, this comes from the fact that the variance 
is dramatically sensitive to large price variations. This indicates that mo- 
ments are not sufficient anymore to capture the informations contained in 
the price and wealth distribution : one must study directly the distributions 
themselves, and more precisely the tails of the distribution P(AW) of the 
wealth variation AW. By the laws of composition of Levy laws, one can 
show that P(AW) decays as ATim+^ f° r large losses (AVF — > — °°)) with 
Wo depending on <p(x, t). We then propose to determine the optimal strategy 
(j)*(x, t) by the condition /^t) = ® ' corres P on ding to a minimization of the 
'catastrophic' risks, since Wo controls the scale of the distribution of losses, 
i.e. their order of magnitude. 

The full derivation of the solution of this minimisation problem will be 
presented elsewhere; here we restrict to the simple 'variational' ansatz where 
4>(x, t) is taken to be a constant 0, independent of x and t, to be optimized. 
Using Eq.(l), we need to calculate the distribution P(AW) of large losses. 
AW can take large negative values when 

• x(t) drops dramatically : the option is not realized (max(x(T)— x c , 0) = 
0) but the bank looses (x(T) — x )(f) due to its hold position. Taken 
alone, this favors — > 0. 

• x(t) increases much above x c : the bank has to produce the share and 
thus loses — (x(T)— x c ) but partially compensates this loss by its holding 
of (x(T) — x )(f), thus resulting in a net loss of —(1 — (f>)(x(T) — x ). 
Taken alone, this situation favors — > 1. 

There is thus a trade-off between the two possibilities, leading to a non-trivial 
optimal 0. 

In order to obtain it, we write the distribution P(AW) of large losses as 
the sum of these two independent processes and get 



wg = (i - v)C-<jf + vc+(i - <f>Y. 



(4) 



where C± describe the 'tails' of the historical distribution: 



P(x,T\x o ,0) 



C± 



— x—xo^ioo 



X — Xo\ l+il 
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V is the probability at t = that the option will be exercised. Minimisation 
of Wo with respect to yields 



[C + P]< 1 r , 

^ [C+P]f+ [C_(l-P)]« lJ 

A more complete treatment with an arbitrary <f>(x, t), in fact leads to the same 
result given by eq.(5), but with a time dependent V = f£° dy'P (y',T\y,t), 
which is the probability at time t at which the stock price is y that the option 
will be exercized (at T). 

If (1 - V)C- < VC+ (resp. (1 - P)C_ > VC+), <\> -> 1 (resp. 0), as is 
natural since then one of the two histories dominate. If both have the same 
weight, the intuitive result <fi* = | is recovered. The resulting optimal 'scale' 
Wo{<jf) is now the correct measure of the risk which must be added to the 
risk neutral price with a coefficient depending on risk aversion. 

The case where /i < 1, for which the mean becomes itself ill-defined, has 
been discussed in . Let us simply mention that the optimal strategy in that 
case is either cff = or 1. Finally, in the border case fi = 2 (asymptotically 
attracted to the 'Gaussian' stable law), with symmetric tails C + = C_, we 
recover (p*(x,t) = V(x,t), which is precisely Black and Scholes's result. 



3 Conclusion 

We have thus proposed a pragmatic and flexible theory of option pricing, 
which naturally generalizes the Black-Scholes results to intrinsically risky 
markets. Our method allows one to decompose in readable way the option 
price into an irreducible part corresponding to a 'fair game' condition plus 
an extra- cost K, which includes the residual risk associated with trading and 
transaction costs. The risk part is minimized in a variational way with respect 
to the strategy. A well defined trading frequency appears, as a trade-off 
between risk and transaction costs. Our precedure is furthermore relatively 
easy to implement numerically, and we hope that it will prove useful to 
professionals. 

We wish to thank J. P. Aguilar, L. Mikheev, J. Miller and Ch. Walter for 
interesting discussions. 
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Figure Captions. 

Fig. I. : Historical distribution of the daily relative variation for the 
MATIF daily quotation variations in the period 1990-1992 and comparison 
with the best fit to a normal distribution (dotted curve). As is by now fa- 
miliar, the empirical distribution has clearly visible 'fat tails' and a sharper 
maximum. For example, the kurtosis of the distribution is equal to 28.1 
instead of 3 for a normal distribution. 

Fig. 2 : a) Comparison between the option price obtained from Eq.(8) 
with A = and assuming independent increments (full line) and from the 
Black and Scholes formula with the historical volatility (dotted line), applied 
on MATIF options with 30 days maturity (with striking price higher than the 
daily quotation x = 100) and put (with striking price lower then x = 100) 
options. As well known, Black and Scholes formula underestimates option 
prices when using the historical volatility. . The difference between the 
two curves is significantly larger than the error bars from our numerical 
implementation of the analytical formula based on historical data. 

b) Inset: 'Implied volatility' obtained by inverting the Black-Scholes for- 
mula. The well known 'volatility smile' is reproduced. 

Fig 3. : a) Risk aversion cost \^R* (dashed line) and total transaction 
cost K (dotted line) associated with rehedging as a function of trading time 
r, for a 30 days option with x c = x . We assume that trading costs 0.05% of 
the 'share' value. The sum of the two curves (full line) represents the total 
extra cost /C that should be considered to adjust the option price when the 
risk aversion coefficient A = 1. Notice that /C has a minimun for r ~ 10 
days (indicated by the cross) which is, for this choice of parameters, the op- 
timal trading time. Total extra cost at r = 10 is .77 which represents an 
appreciable fraction of the 'risk-neutral' price. 

b) Inset: Comparison between the residual risk corresponding to our 
optimal strategy (full line) and the Black and Scholes strategy (dotted line). 
As expected the latter is larger. 
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Inset: Technical sheet. 



In order not to overload the main text with mathematical formulae, we 
give in this separate inset the most important final equations used for our 
numerical implementation, where we have neglected interest rate effects (i.e. 
r = 0). The optimal strategy <p*(x, t) is given by: 

where P(x', t'\x, t) is the probability that the value of x' occurs (within dx') at 

t', knowing that it was x at t < t', (^)( x ,t) >{x',T) is the mean instantaneous 

increment conditioned to the initial condition (x, t) and a final condition 
(x',T), and D(x) is the local volatility. The residual risk TV then reads: 



T 



71* 



+oo 



dxV(x)P(x,t\x o ,0)(f>* 2 (x,t) 



(1.2) 



where TZ C is the "bare" risk which would prevail in the absence of trading 
(0OM) = 0): 



7Z r 



dx(x - x c ) 2 P(x,T\x ,0) 



roo 

/ dx(x — x c )P(x, T\x , 0) 



(1.3) 



Finally, the option price is determined using Eq. (2) as: 



C(x ,x c ,T) 



r+co 

/ dx'P(x',T\x o ,0)(x' - x c ) + \Vn* 

Jx c 



(IA) 



which allows one to obtain numerically the option price, once P(x', t'\x, t) 
is reconstructed, which we do using the histogram of daily variations plus the 
assumption that these daily increments are uncorrelated. 
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